Pseudoscalar glueball and rj — rj' mixing 
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We have performed a dynamical analysis of the mixing in the pseudoscalar channel with the goal 
of understanding the existence and behavior of the pseudoscalar glueball. Our philosophy has not 
been to predict precise values of the glueball mass but to exploit an adequate effective theory to the 
point of breaking and to analyze which kind of mechanisms restore compatibility with data. Our 
study has lead to analytical solutions which allow a clear understanding of the phenomena. The 
outcome of our calculation leads to a large mass glueball AIq > 2000 MeV, to a large glue content 
of the rj' and to mixing angles in agreement with previous numerical studies. 
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I. INTRODUCTION 

Quantum chomodynamics (QCD) is the theory of the 
strong interactions. A particularly good test for under- 
standing its non-perturbative behavior would be to find 
a good description of glueballs and its properties. The 
glueballs are bound states whose valence constituents are 
gluons, the gauge particles of the theory. For this reason 
the glueball spectrum has attracted much attention 
In particular, the existence of the pseudoscalar glueball 
has been a matter of debate since the Mark II experi- 
ment proposed glueball candidates [!]■ It became clear 
later that only one of them the 77(1405) behaves as a wish- 
ful glueball in its production and decays, i.e. it has not 
been produced in 77 , it has comparably large branch- 
ing ratios in J/ip decays and has not been seen to decay 
to 77 [U, H, Besides the 77(1405), other particles be- 
low 2 GeV have been proposed as glueball candidates 
[Hi From the theoretical point of view, while some 
models tend to support this assignment, others, as well 
as, quenched lattice QCD, predict masses over 2 GeV 
[1, 6] 41]. 

For the purposes of this paper we accept the existence 
of at least one pseudoscalar glueball state. Note that 
the pseudoscalar sector is a complex one. On the one 
hand it accommodates the Goldstone nature of the pseu- 
doscalar multiplet, on the other, not totally unrelated, we 
encounter the singlet-octet mixing, which is traditionally 
associated with the resolution of U{1) anomaly. In con- 
stituent models the ideal mixing {6i — tan~^ ^/2) is nat- 
ural, however the 77 and 77' mixing is non ideal. In order 
to describe this phenomenon a complementary picture of 
low-energy QCD, given by an effective Lagrangian where 
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the underlying chiral symmetry is manifest and the res- 
olution of the U(l) anomaly can be implemented, has 
been proposed [sl-flSj. A modification of this effective 
theory can be performed which incorporates the pseu- 
doscalar glueball without loosing the low energy realiza- 
tion of the fundamental properties of QCD fl^ and leads 
to a 77 — 77' — mixing and its consequent phenomenology. 

We proceed here by following this effective Lagrangian 
prescription, but contrary to other authors, we take the 
experimentally known parameters in the meson sector as 
input and leave the glueball parameters, its mass and 
mixing parameters, as unknown. In Section II we re- 
discuss, with modern phenomenology, the 77 — 77' mixing 
in the effective theory approach, to discover that we are 
not able to fit the data. In order to solve the discrepancy, 
in Section III, we incorporate the pseudoscalar glueball, 
following the approach of ref. [T^ . In our approach, con- 
sistency implies that Mq > 1500 MeV. In the next Sec- 
tions we proceed to study the consequences of the theory, 
i.e. J/ip — >■ r]{r]')X, meson radiative decays V — > 
and 77' V-f, and 7^(7;') 27 decays. These calcula- 
tions force us to incorporate phenomenologically addi- 
tional glueball couplings to the octet 77 meson. We are 
able to solve exactly the model with glueball-octet cou- 
pling in the approximation of two mixing angles. Our 
results are compatible with data for glueball masses be- 
tween 2100 < Me < 2300 MeV. 



II. THE CHIRAL LAGRANGIAN 

Guided by symmetry principles, we can build an ef- 
fective Lagrangian describing the low-energy behavior 
of QCD. The relevant degrees of freedom are the Gold- 
stone bosons of the symmetry breaking G ~ SU{3)l <8) 
SU{3)r H = SU{3)v- There are eight pseudoscalar 
Goldstone bosons living in the coset G/H and transform- 
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ing according to 
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In here we assume spontaneously chiral symmetry 
breaking and an implicit integration over the scalar 
mesons. The explicit chiral symmetry breaking is pro- 
vided by a mass term which mimics the one in the QCD 
Lagrangian: 

/:o = ^ (d^U^d'^U) + ^ (MU^ + UM^) . (2) 

We will work with an isospin SU{2) symmetry and the 
mass matrix is Ai — diag(m, m, m^). The eight Gold- 
stone bosons {'K,K,r}%} are collectively represented by a 
non linear parametrization 



U = exp 



' F 



(3) 



with y/2V = P'^Xa or, in term of physical particles 



TT ' 
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The physical (squared) masses are extracted from the 
quadratic term in ([2]) and we obtain 



Ml 



Bfh 

B{m + ms)/2. 



(5) 
(6) 



The parameter B is related to the quark condensate 
through 



{0\qq\Q) = - 



dC. 



QCD 



dr 



dC 
dM 



= -F^B. 



(7) 



and the instantons. As a consequence, the 770 is not a 
Goldstone bosons and its mass should not be given by 
the mass term in i.e. (2Ml + M^)/3. An additional 
mass term should be added. Following the refs [id. |Tl|. 
the U{1)a breaking term involves det J/*-^^ and reads 



F^ a ( det U 

^^-TeivV^'ldiuTt 



3 a 

'2N 



(10) 



where TV is the number of colors, a a dimensionless cou- 
pling and the equation is valid for three flavors. The 
logarithm is essential to avoid the presence of higher or- 
der 770 self-couplings [11]. In ([TU]), we explicitly write the 
N dependence to show that in the large— iV limit, the 
anomaly disappears [Tsi . 

The introduction of the isosinglet rjQ induces a mixing 
with the 778. It is generally assumed that the physical 
particles 77 and rj' should then be a linear combination of 
the two fundamental fields 1431 



cos 9 
sin 6* 



— sin6' 
cos 9 



(11) 



The masses of the two physical states are the eigenvalues 
of the mass matrix for the 7^3 — rjQ system 
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(12) 

In the SU{i)F limit, where all quarks have the same 
mass, i.e. nis — rh, the coupling between the singlet and 
the octet disappears. At this stage, a is an unknown 
parameter. We can eliminate a in terms of the mixing 
angle 9, which can be determined from the two photon 
decays using the formula 



r(77 77) 

r(7r'' — 77) 



M,, 



cos 9- 2V2sin6' 



(13) 



F is the pion decay constant — 132 MeV as it can be 



deduced from the conserved current Af, 
the definition 



-Fd^P" and 



(8) 



With the Lagrangian ©, all Goldstone bosons have the 
same decay constant F = F^^. 

The ninth pseudoscalar boson is not a Goldstone bo- 
son. However, it can be included in a straightforward 
way in the Lagrangian We add to the representation 
the trace with the Tyo meson properly normalized: 



V 



I3/V3. 



(9) 



The matrix U now belongs to U{3). The apparent 
U{1)a symmetry of the effective Lagrangian ^ should 
be broken by an additional term since this symmetry is 
not a symmetry of QCD. The U{1)a symmetry is bro- 
ken in QCD at the quantum level by the axial anomaly 



T{r)' -» 77) ^ 
r(7r" -f-f) 3 \M, 



1 /AV 



sin6l + 2^2 cos 61 



(14) 



The experimental input leads to a first determination of 
the mixing angle, 9 = -20° 

The physical masses are then functions only of the mix- 
ing angle: 



Af2 = 1 
" 3 



4Af|. - Af2 + 2V2{Ml - Ml) tan 9] , (15a) 



1 r 



Af , = 
" 3 



4Af^ - Ml - 2V2{Ml - Ml) cot 9 . (15b) 



Nevertheless, it is not possible to fit the two masses si- 
multaneously [1^. Indeed, as a check, we can eliminate 
the mixing angle 



tan" 9 = 
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Plugging this result, = —11.4°, in (ITSI) does not provide 
the physical masses for the ry and the rj' . Instead, we 
find ihn = 530 MeV and m^/ = 1181 MeV. Indeed, as 
shown by Georgi [ T7| , this mixing scheme cannot provide 
the physical ratio M^,/M^. 

We conclude that with only one parameter (a or 9) it 
is not possible to reproduce simultaneously the masses of 
the r] — r]' system. We investigate in the next Section an 
improvement consisting in incorporating a pseudoscalar 
glueball into the mixing scheme; another, to be discussed 
elsewhere, is the use of different decays constants for non- 
strange and strange mesons. 



III. THE CHIRAL LAGRANGIAN WITH 
PSEUDOSCALAR GLUEBALL 

The motivation for the inclusion of the extra term, 
(|10p , was to implement the axial anomaly in the effective 
Lagrangian. In term of the isosinglet current, the axial 
anomaly reads in the chiral limit 



(17) 



with = (257^755 + S7^75s)/-\/3- Another way to im- 
plement the axial anomaly in the effective Lagrangian 
is to introduce a field Y interpolating the topological 
charge operator G^^G'^" Since we already have a 
pseudoscalar flavor singlet in the Lagrangian, the r^Oi Y 
can then be considered as an auxiliary field introduced 
via the term 



Ca 



N \ \dciW 



1, 



-Y' 



(18) 



which is equivalent, as can be shown using the equations 
of motion, to (|10l) . But the operator G^^G^*^ may also 
interpolate a pseudoscalar glueball [31 • We split the Y 
field into an auxiliary field, r^aux, describing the Tyo and 
another, g describing the glueball. We add a kinetic term 
and, for the sake of completeness, a mass term associated 
also to the pseudoscalar glueball [l3] 



Ca =«(?/aux + g) { In 



detU 
det C/t 



2 ''l^aux 



(19) 



The first term induces a coupling between 770 and the 
pseudoscalar glueball but we do not have any coupling 
between the glueball and 773 • The mass matrix has then 
the simple form [3| 



Ml,, 



AM^ - 
-2^/2{Ml 




Mi 
-Ml) 



-2V2(Af2 




Where we have defined a = 48/ci, /3 = (4/i^)^6/(c3Ci) 
and 7 = C2/C3. The eigenvalues of the matrix represent 



the mass of three physical states, 77, t]' and a third pseu- 
doscalar state Q. In order to simplify the relations, we 
introduce the following notation 



'W 



0' 



/? 7; 



with 



W ^ - (AM 



K 



Mi 



Z^-^-^{Ml-Ml), 
y = i {2Ml + Ml) . 



(21) 

(22) 

(23) 
(24) 



The mass matrix M/^^g is diagonalized using a rotation 
matrix R 



(25) 



TW^ = diag(M2,M2,,Af2) ^j^j^ ^.^^^ ^Yie unknown mass 
of the third, mainly gluonic, state. The rotation matrix 
R collects the eigenvectors of the transformation between 
the pure states and the physical states. In ref. fis'l, the 
authors used the eigenvectors R to diagonalize a matrix 
linear in the masses but the chiral Lagrangian leads to 
relations quadratic in the masses. Obviously, there is no 
contradiction between these approaches. 

The knowledge of R determines the decay properties of 
the physical states. The matrix relation provides 6 
independent relations since the matrix is symmetric. We 
have three unknown parameters in A^sog (Q!,/3,7) and 
one in Ai, the mass of the third pseudoscalar particle 
Mq. If we could find a rotation matrix in terms of two 
mixing angles, all these quantities could be determined. 
This hypothesis is often considered in the literature where 
the rotation matrix is parametrized with two angles, one 
for the rotation between 770 and Gluonium and a second 
angle for the rotation between r/o and 773 [T9| . However, if 
we accept the existence of a real gluebal state, i.e. 7 7^ 0, 
it is not possible to obtain a matrix of the form of ((2T|) 
for Aisog with only two angles. As will be shown in Sec- 
tion |Vl the two mixing angle scheme is recovered in our 
description if we incorporate an octet-glueball coupling. 

Without any assumptions on i?, i.e. with the more 
general three angle Ansatz, we can only determine the 
parameters a,/3,7 as functions of Mq. For this purpose, 
we have to equal the three rotation invariants given by 
the coefficients of the characteristic polynomial, P{X) = 
X^—TX'^+SX—D, of the matrices. The three invariants 
for are 



D = M^M^,Ml, 
f2 n j-2 



S = M^M^, 



T 



Mi. + Mi 



Ml,, Ml 
+ Ml. 



MlMl 



(26) 
(27) 
(28) 



Those three quantities are function of Mq since we take 
the physical masses for the rj and 77'. 
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It is now easy to extract the values of the parameters 
in terms of the known quantities 



7 



= W 



1 

Z2 



(M^3 _ ^SW -D) (29a) 



(7 + T4^)(T-T4^) - (^2 + 5 + 7^), (29b) 
r-(VK + y + 7). (29c) 



Only if > our system will have a solution. This con- 
dition restricts the allowed values for the glueball mass 
Afe- The equation /J^ = is quadratic in Mq leading to 
two solutions given by 



and X G {8,0,g}. Sometimes we use the strange and 
non-strange components of the eigenvectors which are ex- 
pressed in term of the previous comp onents by a rotation 
with the ideal angle (cos 9i = -^/l/S) 



cos ti; 

- sin ( 




(31) 



We use the following convention to denote the compo- 
nents in the non-strange— strange basis: 



\P) = Xp\ijq) + Yp\i^,) + Zp\G), 



(32) 



= W 



(30a) 
(30b) 



which are shown as a function of Mq in Fig. [T] The 
bounds on Mq are the extension of Georgi's bound 
in the case of a third pseudoscalar particle. 




FIG. 1: as a function of Me- 

In ref. the r] mass was adjusted, since the lower 
bound of does not depend on the t]' mass, to have 
a positive for Mg ~ 1400 MeV. Their aim was to 
accommodate the 77(1405) in the theory (43|. In our study, 
we leave the Mq as a parameter and therefore the > 
condition implies Mq > 1500 MeV. 

IV. J/ip DECAYS 

The theory we have just described contains an unique 
parameter, the glueball mass Mq, out of which we can 
extract many consequences which are observable. We 
will center our attention in the J/ip to 7], rj' , and 7], rj' 
to two photon decays. These decays are described in 
terms of the components of the eigenvectors, the rows of 
R in (j25p and correspond to the mixing parameters for 
the physical states They will be labelled as and are 
defined through \P) = Y.x^x\^^) ^i^h P e {?7,7]',e} 



where Vg = Z. 

The data that we attempt to describe have been taken 
from the Particle Data Group compilation [20| . 

The radiative decays of the J/tp into 77(77') take place 
through the anomaly (0|G^,yG^'^ 177(77')) and their branch- 
ing ratio is given by 



r( j/y, ^ rj'j) 
r( j/^ ^ 777) 



M2 



4.81 ±0.77 



(33) 



Other J/ip decays which may probe the strange and 
non-strange quark contents of the 77 and 77' are the 
ones producing (f>, and oj{p) respectively. The processes 
J/ip ^ ^{11') P violate G parity and isospin. They pro- 
ceed through the exchange a of virtual photon 2l| and 
we find 



T{J/i> ^ 77'p) 

r( j/t/- ^ 77p) 




0.54 ±0.16 



(34) 



The pseudoscalar meson momentum in the center-of- 
mass is k"^ = A(M2, ,M|,,M^)/(2Mj/^) defined in 



terms of 



\{x, y,z) — \/a;2 + 7/2 + — 2xy — 2yz — 2zx. (35) 

A standard approximation (^//^ ^ MpMy) for the 

2 



momentum is kp ~ Af 



-j/^(l-(Af2+M2)/Mj/^)/2[21|. 

The processes J ^ v{v')<^ ^ind J/ip vi''l')4' pro- 
ceed again through the exchange a of virtual photon but 
also via OZI processes. In QCD, three gluons are emitted 
from the J/ijj and give rise to two light quark- ant iquark 
pairs. In our effective approach, the degrees of freedom 
are the mesons and the interaction is modelled by the La- 
grangian eafSf,ud°'d^{Tt'{V''r)) where T'^ is the J/t/j field 
and all the light axial- vector mesons are collected in the 
matrix V in analogy with the representation in Eq.Q. 
Both the isospin violating and OZI processes have to be 
taken into account for a complete description. Since we 
do not want to add at this stage one more parameter 
we compare branching ratios. Assuming an ideal mixing 
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between w and 
butions 



we find for the sum of the two contri- 



T{J/iP -> r]uj) 



r(j/^ ^ 770) 



0.105 ±0.024 




0.53 ±0.15 



(36) 



(37) 



With the momenta kp defined as above. Within this 
model, the formulas for the p and oj decays are similar. 
Since Mp ~ M^j, the phase space is almost equal and the 
model predicts the same value for their branching ratios. 
However, the experimental data give a factor 5 differ- 
ence. This discrepancy is eliminated if one incorporates 
the contribution of more sophisticated decay processes 
like the double OZI processes [H,!!!]- 

We display these ratios together with the experimental 
data (in gray) in Fig. dH and El 




FIG. 2: T{J/i) ?7'7)/F(J/i/' ^ V't) as a function of Me in 
the scheme without octet-glueball coupling (Sect. IV). 




FIG. 3; r(J/V' ri'p)/r{J/i) rjp) as a function of Me in 
the scheme without octet-glueball coupling (Sect. IV). 




FIG. 4: T{J/il} rj'u)/V{J/4j rfu)) as a function of Me in 
the scheme without octet-glueball coupling (Sect. IV). 




FIG. 5: V[J/ip ri'cj))/r{J/ip -qcj)) as a function of Me in 
the scheme without octet-glueball coupling (Sect. IV). 



In order to show the amount of mixing in a specific case 
we give the mixing matrix for Me — 2000 MeV (which 
corresponds to the red diamonds in the figures): 



0.9874 
0.1492 
-0.0534 



0.1107 
-0.4085 
0.9060 



-0.1133 
0.9005 
0.4198 




(38) 



It is worth mentioning that, if there is no coupling 
between the glueball and the octet [see ea. (l20t ]. and if 
we force the t] and rj' masses to their physical values, the 
particle with the most gluonic content is the rj' ! This 
statement, which remains valid for any Mq, is clearly in 
contradiction with the usual assignment for the 77', since 
the most gluonic particle should be Q. This is a clear 
indication of the need for improvement. Moreover, as 
shown in Fig. [5][S1 this mixing scheme cannot explain the 
data on the J/V' decays except for J/V' — > ri{ri')uj. 

The model thus far developed does not capture the 
physics of the pseudoscalar sector with our philosophy 
consisting in fixing the meson masses to their experimen- 
tal values [l^l- One mechanism to cure these incompati- 
bilities is to introduce a coupling between the octet meson 
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and the glueball. 



V. COUPLING BETWEEN OCTET AND 
GLUEBALL 

The fact that the strange quark mass is heavier than 
the up and down quark masses leads in model calcula- 
tions [2^ [IJ] which implement QCD in a confined sce- 
nario to a non vanishing coupling between yyg and G. In 
Fig. E] we show a possible diagram which contributes to 
the mixing. 



Ref. 



decays 



n 



W30D0D0000000' 



[27] P(V) ^ y(P)7 41.5 ±1.2 -13.2 ± 1.2 0.04 ± 0.09 

[28] P{V)-^V{P)'Y 41.3 ±0.7 -13.4 ±0.7 0.04 ± 0.04 

[28] J/ijj-fVP 45 ±4 -13.4 ±0.7 0.04 ± 0.04 

[19| J/ip^VP 44.5 ±4 -9.7 ±4 0.28 ± 0.21 

[29| P(T/) ^ y(P)7 40.4 ±0.6 -14.3 ± 0.6 0.12 ±0.04 

TABLE L Summary of the recent work on the »7 — 77' — {glue) 
mixing. The second column is given hy 9 = (p — 6i. 



The particular Ansazt gives the same mixing scheme 
for the 77 as in (fTTj) . Hence the relation (l42aD is equivalent 
to the previously derived relation (|15al) . This theoretical 
estimate 9 = —6.4° is Me— independent and agrees to 
Icr with the recent numerical study of Escribano |T9|. 
e = (-10.2 ±4.3)°. 

The second of the equations (|42)l leads to the same 
constraint as previously found, 



FIG. 6: Mechanism within QCD of Octet-Glueball mixing. 

We here take a phenomenological point of view which 
consists in enlarging the mass matrix to incorporate this 
coupling by means of a new parameter (5[45|. 



'W Z 
Mlog= I Z Y + a P 



(39) 



M^>W 



W 



(43) 



In the presence of an octet-glueball coupling we do not 
find the unphysical branch {Mq < M^) but only the 
bound Afe > 1-5 GeV. The second angle is displayed in 
Fig. [7] as a function of AIq. It is worth mentioning that 
we expect a lower mass bound when including different 
decay constants for the mesons. 



In order to reduce the number of unknowns, here 
one more since we added S, we choose next a rotation 
parametrized with only two angles, 



R = 



with the convention 




cos 9 — sin 6* ' 
sin 9 cos 9 
1, 



(40) 









0.2 


i) 




(41) 


0.1 












This particular form for the eigenvectors assumes no glue 
content in the ry wave function. The discovery of a p he- 
nomenological evidence for glue content in the ij' [26| led 
to the introduction of additional angles in the rj — rj' mix- 
ing scheme ([TT|) . Generally only a second angle is added 
and the scheme ([Ml is assumed [19, 27 29]. We present, 
in Table HI the summary of the most recent studies on 
this topic. 

From the six equations of the matrix relation A^log = 
R^M^R, we find 



tan^ 



cos 



W 



M2 

77 



Z 



w ^ 



Ml 



(42a) 
(42b) 



FIG. 7: Me-dependence of sin 



It is also possible to extract analytically the values of 
the four couplings in terms of the physical masses 



Y 



7 = Ml,+Ml-W 



a = T-W 



M2' 



6' 



^'K'-7)(7-M|) 



Z2 



{W - M2)2 



(IF- 



Z2 



M2)2 



(44) 
(45) 
(46) 

(47) 
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We see that 7 G Wl^M^]. 

In Fig. m ini and [tni we compare the calculation with 
the data. The experimental values for the ratios are dis- 
played in gray. 

As explained in the previous section, our formulas for 
p and LJ decays are similar and therefore not consistent 
with the data. We are thus not able to fit simultaneously 
the p and uj decays of J/ip. However, in the case at hand, 
with an octet-glueball coupling, the (j) decay of the J/ip 
is consistent with the p decay within the interval 



2.1 GeV < Mq < 2.3 GeV 



(48) 



In terms of the glue content of the rj' , the interval pS)) 
reads 0.38 > Z^j' > 0.30. The description of the w decay 
would require a lower glueball mass. 

There is no radiative decay J/tp — ^ 777 since, in the 
two angle rotation scheme, there is no glue content in 
the 77 and therefore the corresponding branching ratio to 
the rj' cannot be defined. In the future one might want 
to study a three angle rotation scheme which however 
requires numerical treatment. 

In order to see the amount of mixing in a specific case 
we show the mixing matrix for Mq = 2200 MeV (repre- 
sented by red diamonds in the figures) which corresponds 
to (j) = 35.7° (The other angle 6 = —6.4° does not depend 
on Mq): 



ri' 



0.9938 0.1114 
-0.0904 0.8065 0.5842 
0.0651 -0.5806 0.8116 




(49) 



2.2 2.4 

M„ (GeV) 



FIG. 9: r(J/'i/' — > r]'uj)/r{J/^l! tjuj) in the two angle scheme 
(Sect. V). 




0.7 - 

!=■ 0.6 - 




2.2 2.4 

M„ (GeV) 



FIG. 8: T(J/ip -> ri'p)/r{J/tp — > -qp) in the two angle scheme 
(Sect. V). 



FIG. 10: r(J/V' 
scheme (Sect. V). 



r]' (f)) /V{J — !> rj(f)) in the two angle 



use the following relations and we quote the experimental 
values: 



r(w 777) 
r(r?^ ^ Pi) 

T{p rji) 

r(0 ^ v'l) 



= 3 



X 



n' 



X 



1.58 ±0.43 

1.35 ±0.24 

2 



M2 



M2 



M2 



M2 



Y, 



Ml - A4 
Ml - M2 



(4.78 ±0.25) 10" 



(50) 
3 

(51) 
(52) 



VI. MESON RADIATIVE DECAYS 

In this section, we explore the meson radiative decays 
V — > Tlii]')"/ Slid v' ~^ Vi- The interacting Lagrangian 
modclhng those decays is ea,3^^i^"^9^(V''P) with F"^ 
the field strength for the photon. In particular, we will 



We display in Fig. [TTl [TH and [T3] the radiative decays 
between pseudoscalar and axial- vector involving cj, p and 
(/) decays respectively in both schemes. 

In view of the results, we can discard safely the first 
model without octet-glueball coupling. When such a cou- 
pling is introduced, the data for p and a; are consistent. 
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FIG. 11: Radiative decays r(?;' — >• u]'y)/Y'{uj — >■ 777) in both 
schemes (dashed line: without octet-coupling and solid line: 
with octet-glueball coupling). 



FIG. 13: Radiative decays V{<f) -> 'q''y)/V((f) 777) in both 
schemes (dashed line: without octet-coupling and solid line: 
with octet-glueball coupling). 



- without coupling 

- octel-G coupling 




2.2 2.4 2.6 2.8 
M„ (GeV) 



FIG. 12: Radiative decays V(ri' — )■ p'y)/r{p — >• 777) in both 
schemes (dashed line: without octet-coupling and solid line: 
with octet-glueball coupling). 



Moreover, the allowed range for physical glueball mass 
Mq lies in the same range as for the J/ip strong de- 
cays dm). 

The (j) radiative decay does not fit the data. We see no 
mechanism to lowest order to cure this problem. 



VII. DECAYS INTO PHOTONS 

The strong and radiative decays allowed us to discard 
one model and forced us to introduced an octet-glueball 
coupling. In this section, we now study the decays into 
photons. Electromagnetic decays are more sensitive to 
the decay constants. In our model we use the same decay 
constant for all pseudoscalar particles, nevertheless we 
expect to have good quantitative results. 

In order to calculate the decays of the pseudoscalars 
into two photons, we add the Wess-Zumino-Witten 



(WZW) term 



C 



WZW 



(53) 



Q2 _ diag(4/9, 1/9, 1/9) is the matrix of the quark 
squared charges and F^^y is the field strength for the 
photon [4^. This term only couples quarks to photons 
since the gluon does not carry electric charge. We obtain 
for the branching ratios. 





> 77) 


\ 


( 


r(7r0 - 


-> 77) 






r(V- 


-> 77) 


\ 


( M^, 


r(7r0 - 


77) 






r(e - 


-> 77) 


-\ 


( 


r(7rO - 


77) 







(54a) 



(54b) 



. (54c) 



All the three branching ratios in Eqs. (IM)) . can be re- 
cast in the form 



r(p ^ 77) 



Ah 



4. 



r(7rO 77) \Aho 
The experimental values for these coefficients are [l 



(55) 



Cri = 0.944 ±0.040, 
c„. = 1.242 ±0.027. 



(56) 
(57) 



The 77, rj' decays into two photons are shown in Fig. [HI 
and 1 151 for the two Ansatze used in the previous sections. 
We notice that it is not possible to reproduce the data 
without glueball-octet coupling. In ref. lj|, the authors 



used the value of the 77 mass as a parameter to accom- 
modate their model to the data. If we use the physi- 
cal mass of the 77, the branching ratio for the 77 is quite 
Me— independent and not in agreement with the data. 
Tlie 77' decays is neither in agreement. 

Our second parametrization, the two angles scheme 
with octet-glueball coupling, leads to an 77 with no glue 
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1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 



FIG. 15: (blue) and c^/ (red) in the two angle scheme. 



cate, that the inclusion of different decay constants for 
the pseudoscalar meson octet, might change considerably 
the 6 mass range, leading to lower allowed mass values. 
However, we can conclude safely that Mq > M^i. 

Having said this, we can present our model predic- 
tions for decays involving this third pseudoscalar and 
its decays. In Fig. [1^] and [T7] , we plot as function of 
Me, the branching ratios r{J/ip — ^ Qp)/T{J/ip — )■ rjp) 
and r(J/i/> ^ ej)/r{J/ij -> 7/7). in Fig. M 
r(0 — 77)/r(7r° — > 77). We see that these observables 
are non overlapping. The J/ip branching into X — p or jj 
in the @ mass range are very small while X — j branching 
ratio is large. Unhappily there are no data in this mass 
range. 

0.7 1 > > > > > > > > 1 



0.6 




1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 

Mg (GeV) 



FIG. 16: r{J/ip Qp)/V{.J/iP rjp) as a function of Me. 



content and therefore the value of its branching ratio re- 
mains the same, away from the data. However, the rj' , 
gets a large glue content, leading to a branching ratio 
within the data in the allowed range for Mq, Eq. 
Since these electromagnetic interactions strongly depend 
on the decay constants, we hope to improve the rj de- 
cay into photons by using different decay constants for 
the octet [iTj. In summary the results based on electro- 
magnetic decays strengthen our conclusions based on the 
analysis of strong decays. 



VIII. DECAYS INVOLVING 9 AND 6 DECAYS 

The above discussion has fixed not only our theoret- 
ical scheme but also our parameters. We aim now at 
predictivity. However, we must keep in mind that our 
calculation is a first order calculation (equal decay con- 
stants for all the pseudoscalar meson octet) and therefore 
we expect changes at higher order. The present results 
seem to indicate however, that we are obtaining a satis- 
factory mixing scheme but that we should not trust our 
Q mass range quantitatively. Primitive estimates indi- 




1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 

M,;, (GeV) 



FIG. 17: r(J/7/> Q'y)/r{J/i! r)'-)) as a function of Mq. 



IX. CONCLUSION 

We have performed a dynamical analysis of the mix- 
ing in the pseudoscalar channel with the goal of under- 
standing the existence and behavior of the pseudoscalar 
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0.1 




FIG. 18: ce as a function of Mq. 



glueball. Our philosophy has not been to aim at pre- 
cise values of the glueball mass but to exploit an ade- 
quate effective theory to the point of breaking and to 
analyze which kind of mechanisms restore compatibility 
with data. Our study has lead to analytical solutions 
which allow a clear understanding of the phenomena. 

Let us summarize the main findings of the present in- 
vestigation. Starting from an effective Lagrangian for- 
malism, which incorporates the pseudoscalar glueball, we 
try understand the rj — rj' mixing phenomenology and the 
dynamics it implies. Our approach differs from others in 
the same line [l3, 1131 it that it takes the meson masses 
from experiment and only leaves the glueball parameters 
to be determined. Moreover, we do not proceed by fitting 
but find analytical solutions to the mixing problem. 

The calculation of the J/ if) decays in the initial ef- 
fective Lagrangian is unsuccessful in the explanation of 
the data. Implementing, in a phenomenological way, the 
octet-glueball couphng inspired by QCD, leads to an ex- 
act solution in terms of two angles which fits the data 
for large glueball masses Mq > 2000 MeV and leads to 
a vanishing glueball component of the 77 and a large one 
for the T]' . Our results are compatible with the mixing 
schemes of KLOE and Escribano [l^, [sj] which reinter- 
preted in terms of our mass matrix lead to octet-glueball 
coupling. The chiral Lagrangian to first order in ex- 
tended to include the glueball predicts a mixing angle 9 
compatible at la with previous numerical studies and a 
Me— dependant angle (f) also compatible for a wide range 
of the pseudoscalar gluebaU mass. Our study is a strong 
theoretical justification of the previous analysis of the 
r] — 77'— glue system. 

The 27 decays teach us that the WZW photon cou- 
pling is sufficient to explain the data provided that we 



incorporate an octet-glueball coupling in the model. This 
supports our conclusion based on strong decays. 

Our analysis therefore leads to a new dynamical sce- 
nario which needs to be constructed from the point of 
view of an effective Lagrangian theory. Within this 
scheme we have obtained a compatibility with the data 
for large gluebah masses 2100 MeV < Me < 2300 MeV, 
and large glueball component for the 77'. This large 
glueball mass raises the question of the inclusion in the 
mixing scheme of higher resonances. Indeed, the pseu- 
doscalar spectrum is rich of resonances around 1-2 GeV. 
In this work, we only considered a third gluonic state in 
addition to the usual r]g, and tjq but at this high energies, 
it could be relevant to include other fields in the mix- 
ing scheme such as multiquarks states (36l - [39l |. However, 
even if our result clearly indicates a large glueball mass, 
this has to be taken with a pinch of salt since we have 
used in our scheme the meson couplings as = Fk- If 
we naively relax this assumption followiiig the methods of 
the current algebra schemes [2l|, Isil - lssl . [35l | we can show 
that the lower mass limit decreases considerably and that 
we can expect Mq < 2000 MeV. Moreover, as explained 
in Sec. II, the Chiral Lagrangian at leading order leave 
a lot of room for improvement. We chose to improve it 
with a glueball field but we learn from [s^, that the 
room for the glueball, and hence for other multiquarks 
configuration, is drastically reduced when going at next 
to leading order. 

Our analysis leads to consequences of for further stud- 
ies. We need to construct the effective Lagrangian that 
incorporates octet-glueball coupling. Moreover, we have 
to describe in the effective Lagrangian approach the 
Ftt 7^ Fk dynamics. Certainly our analytical solutions 
are in some aspects naive, but certainly they allow a clear 
interpretation of the phenomena and may serve to test 
all these improvements. 
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It is believed though that unquenched calculations will 
lower this mass [H, Bl- 

This parametrization is an oversimplification used for the 

purpose of illustration [3 • 

The 7^(1405) was at that time the t(1440) . 

The solution to the mass equation has a second branch for 

< A/e < 500 MeV. This mathematical solution is not 

compatible with the data and therefore has no physical 

reality. 

There are many Lagrangian terms which would produce 
this coupling, however at present we do not see any ar- 
gument to choose one. 

In this case, we make a difference between u and d quark 
since they carry a different electric charge. 
One must also keep in mind that gluons can couple to 
photons through quark boxes, althoug we expect this 
mechanism to be less effective then the SU(3) breaking 
in the currents. 

For the present Q mass range the decay into Q — <j)is not 
allowed. 



